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o | Abstract 

CD 

ry^ The bending energy of any freely deformable closed surface is quadratic in its curvature. In the 

absence of constraints, it will be minimized when the surface adopts the form of a round sphere. If 
the surface is confined within a hypersurface of smaller size, however, this spherical state becomes 
inaccessible. A framework is introduced to describe the equilibrium states of the confined surface. It 
is applied to a two-dimensional surface confined within a three-dimensional hypersphere of smaller 
radius. If the excess surface area is small, the equilibrium states are represented by harmonic 
deformations of a two-sphere: the ground state is described by a quadrupole; all higher multipoles 
03 are shown to be unstable. 

B 

^ l 

> 

^ Twenty five years ago, Nickerson and Manning examined how an elastic curve minimizes its energy when 

OO its freedom to bend in space is constrained by a two-dimensional obstacle [HE]- Their motivation was 

^] to understand the wrapping of DNA, a semi-flexible polymer, around histone cylinders but the question 
is also relevant to an understanding of the confinement of DNA within viral capsids (see, for example, 

S 0) 

j> Whereas a space curve minimizes the Euclidean distance between two points by negotiating obstacles 
along surface geodesies, an elastic curve generally will not. This is because the relevant energy is not 
arc-length but the bending energy of the curve in three-dimensional space, quadratic in the spatial Frenet 
curvature k [4J. This energy is always greater than its surface counterpart, quadratic in the geodesic 
curvature n g , characterising how the curve bends within the surface. By projecting the curvature vector 
of the space curve along the surface and normal to it, its bending energy can be decomposed as a sum 
of two contributions: k 2 = n 2 g + [5]; whereas the intrinsic K g is insensitive to the three-dimensional 
environment, the normal curvature K n will register how the surface itself bends in space. When the 
surface is a sphere, K n is a constant and thus the spatial bending energy will be minimized by surface 
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geodesies. The snag, however, is that the boundary conditions, and in particular the closure of the curve, 
tend to be incompatible with geodesic behavior. 

In |F] the equilibrium states of a closed elastic loop confined within a sphere were examined. It was 
found that, if the loop is short, there is a unique ground state; it is completely attached and it exhibits 
a two-fold dihedral symmetry, shown in Fig. [T| This state is stable; all of its excited counterparts are 
unstable. The bending energy was also found to depend-in a surprisingly sensitive way-on the length of 
the loop. Contradicting naive expectations, neither the energy nor the forces transmitted to the sphere 
increase monotonically with loop length. Local maxima are associated with the incommensurability of 
loop length with the equatorial circumference. 




Figure 1: Equilibrium state of an elastic loop with two-fold symmetry 

In this paper, higher dimensional analogues of a confined elastic loop will be examined. In particular, 
the focus will be on confined two-dimensional surfaces. The surface energy is now the two-dimensional 
bending energy, quadratic in its curvature (see, for example, [3 [HJ QUI [HI D-2J). In the absence of 
constraints, it is simple to demonstrate that the closed surface minimizing this energy will be a round 
sphere. 

The obvious higher-dimensional analogue of a confined loop is a two-dimensional closed surface confined 
within a sphere of smaller area. This is a problem of considerable interest from a biological perspective- 
where the morphology of cellular compartments may involves the confinement of one membrane within 
another-and it has been explored recently by Kahraman et al. (T3J [T4] . In distinction to the loop, which 
attaches everywhere to the sphere, the confined surface must detach somewhere, folding into the cavity. 
The interesting behavior occurs within these interior folds and along the boundaries of the regions of 
contact, not within the regions of contact themselves. 

A more appropriate analogue of the confined loop is provided by locating the surface in a higher di- 
mensional environment. The simplest possibility is to embed the two-dimensional surface in a four- 
dimensional Euclidean geometry, confining it within a three-dimensional hypersurface, and, in particular, 
a hypersphere of smaller radius. While this may appear a somewhat recondite exercise from the practical- 
application-driven-point of view of condensed matter, geometrically it is the natural generalization of 
the confined loop. Mathematically, it also can be viewed as a natural generalization of the Willmore 
problem involving competing length scalesjTl [S]. Indeed, when the constraint on the area is relaxed, it 
reduces to the Willmore problem as formulated on a three-sphere, S 3 . 
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To address this problem, the constraint is imposed by introducing a local Lagrange multiplier. This 
generalizes the variational framework developed in [6J. The advantage of this approach is that the 
multiplier quantifies directly the loss of Euclidean invariance associated with the constraint; it becomes 
a source for the stress induced in the surface, its value at any point identified with the local normal force 
associated with the constraint. 

If a surface of area greater than Air is confined within a unit hypersphere it will generally buckle into a 
non-spherical shape. Unlike its loop counterpart, however, this problem is not tractable analytically. In 
this paper, it will be treated perturbatively in the regime where the surface area is marginally greater 
than An. In analogy to a circular loop confined within a sphere, it is found that the confined equilibrium 
shapes are quantized. These states are represented by spherical harmonic deformations of a two-sphere. 
The ground state is a quadrupole, exhibiting a three-fold degeneracy; all higher multipoles turn out to 
be unstable. 

The paper is organized as follows: In section [2j the variational framework is introduced. The Euler- 
Lagrange equations for the surface are derived and the constraining force identified. Whereas the local 
confining force typically is not constant, it will be if the confinement is within a hypersphere. In section 
[3] the weak confinement of a two dimensional sphere by a spherical hypersurface is examined. The 
paper concludes with a brief discussion and suggestions for future work in section |4j A number of useful 
definitions, identities and derivations are collected in a set of appendices. 



2 Confining surfaces by hypersurfaces 

Consider a surface Y N of fixed area constrained to lie on a hypersurface T, N+1 in an iV + 2-dimensional 
Euclidean space [15jQ The surface and hypersurface embeddings in Euclidean space are described 
parametrically by the maps Tjy '■ (s 1 , • • ■ ,s N ) — > Y(s 1 , • • • ,s N ) and £jv+i : (u 1 , • • • ,u N+1 ) — > 
X(m 1 ,-- - ,u n+1 ) respectively. The embedding of the surface on the hypersurface is described by a 
third mapping T NiN+1 : (s\ • • • ,s N ) ->■ (U 1 ^ 1 , ■ ■ ■ ,s N ),--- ,U N+1 (s\--- ,/)). The three map- 
pings are related by composition, Y = X o U. While the problem of specific interest will be in a 
two-dimensional surface embedded in a three-dimensional hypersphere, the dependence of the confine- 
ment process on the dimension as well as the confining geometry is not without interest and involves no 
significant additional calculational cost. 

If the surface is freely deformable its bending energy is given by [TU] 

H=\J dA N (V 2 Y)\ (1) 

where V 2 is the Laplacian on Tn and An its area, both of which are constructed using the induced 
metric on IV, g a b, as described in the appendix]^] 

In order to enforce the condition that lies on Sjv+i, a term enforcing this constraint is added to the 

1 Points here are represented x = a^x,, with Xj the orthonormal basis in R w+2 , (i = 1, . . . , N + 2). 

2 lt is more usual to express H in terms of curvatures. This is done explicitly in Appendix A For surfaces of dimension 



higher than two, H is not the most general bending energy of a closed surface. For example, if TV > 3, there is an additional 
"Einstein-Hilbert" bending energy, proportional to the integrated scalar curvature. While this energy is topological when 
TV = 2 it is not for higher dimensional surfaces. 
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energy H[Y] given by Eq. ([I]): 



H C [Y, U] = H[Y] + J dA N \{s) ■ \Y{s) - X(tf(s))] + p. (^j dA N - A^j 



(2) 



which involves a vector-valued Lagrange multiplier A defined along on T^. This constraint will break 
the manifest translational invariance of the energy H. A second term enforces the constraint fixing the 
area. It involves a Lagrange multiplier p which is constant. 

The response of H c to a change Y — > Y + 5Y in the embedding functions can be cast in the form 

5 Y H C = [ dA N (V a P + A) • 5Y , (3) 



where V a is the covariant derivative on compatible with the induced metric g ab and the surface 
stress f a is given by 

r = i a B -pg ab e b ; (4) 
it involves a contribution associated with the bending energy given by [16] 



f a B = Kj \K ahI - ^g ab K T ^ e b - V a Kj n 1 , (5) 

as well as a tension associated with the constraint on the area. The simplest derivation of Eq. ([3]) 
involves a straightforward generalization of the method of auxiliary variables introduced in [17]. The 



details are provided in Appendix B Here e a = dY/ds a , a — 1, • • • , N are the tangent vectors to the 
surface adapted to the parametrization; the corresponding induced metric is given by g ab = e a ■ e b . The 
vectors n/, I = 1,2 are any two mutually orthogonal unit vectors into IR^ 2 normal to and, for 
each value of I, K^ 1 = e a ■ dijx 1 is the extrinsic curvature (a symmetric surface tensor) associated with 
the rotation of the normal n 1 towards the surface as parameter curves are followed along the surface. 
K 1 = g^K^ 1 is its tracej^] Surface indexes are lowered and raised with the metric tensor and its inverse, 
g ab ; normal indexes are lowered and raised with the Kronecker delta. The covariant derivative defined by 
V a K ! = (5 1 jV 'a + uij j)K J , involves the normal connection cu a IJ defined by u a IJ = n 1 -d a n J = —ui a JI , 
associated with the freedom to rotate the normals among themselves. Thus the analogue of the Gauss 
equations for the normals is dan 1 = K a bI e b — co a I J ii J . This description is completely invariant with 
respect to local rotations of the normal vectors. Confinement will select the normal to the hypersurface 
as one of these normals, thus breaking this invariance. 

In equilibrium, Eq. ([3]) implies that 

V a P = -A. (6) 

Thus, in the presence of the constraint, the stress in the surface is not conserved; the Lagrange multiplier 
A is identified as the external force associated with the constraint [§]. 

The corresponding variation of H c with respect to U A (s) is given by 

5uH c = - [ dA N X-E A 6U A , (7) 



3 Note the identity for the energy density (V 2 Y) 2 = K^i. 
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where = d~K/du A , A — 1, • • • , N + l are the N+ 1 tangent vectors to the hypersurface adapted 
to the parametrization defined by the coordinates u A . In equilibrium, A ■ = 0, and thus the force on 
the curve associated with the contact constraint always acts orthogonally to the surface. One thus can 
write A = An, where n is the unit vector normal to the hypersurface, £jv+i, and thus 



-An 



(8) 



Using Eq. (|4]), a straightforward calculation implies that the divergence of f a is normal to the surface 
and may be decomposed accordingly: 

V a f = eX, (9) 

where the Euler-Lagrange derivatives ej of the bending energy along the normal directions n J , are given 
by 

e I = - V 2 K! + ^{KjKj - 2K Iab Kf)K J + yXj , (10) 

where V 2 = g ab V ' a V&. The two tangential Euler-Lagrange derivatives V a f a • ej vanish identically, a 
consequence of the fact that the Hamiltonian depends only on the surface geometry so that tangent 
deformations get identified with reparametrizations. 

The surface normals can always be rotated so that one of the two coincides with the hypersurface 
normal n. The second normal (1) is then uniquely fixed (modulo a sign) by the conditions, 1 ■ n = 0, 
and 1 • e a = 00 With the choice n x = n and n 2 = 1 it is possible to write 



V Q f a = e n n + e x 1 . 

Comparing Eq.(12) with Eq. ([8]), the Euler-Lagrange equation 

£l = Q 

is identified. The corresponding projection onto n determines the normal confining force A: 

A = -e„ ■ 



(12) 



(13) 



(14) 



Both the Euler-Lagrange equation (13) and the confining force (14) are constructed out of the two un- 
constrained Euler-Lagrange derivatives. While the latter is completely determined by the local geometry, 
what this local geometry is will, of course, depend on the global behavior of the confined surface. 



Eqs. (13) and (14) are still not very useful in their present form. To facilitate their interpretation as 
well as their application, it is useful to cast them explicitly with respect to variables adapted to the 
environment of the fixed confining hypersurface. 

Let n ab represent the extrinsic curvature of T NiN+1 embedded in T, N+1 (associated with the rotation 
of the normal 1), and Kab the extrinsic curvature of embedded in M. N+2 (associated with the 



rotation of its unique normal n) as described in Appendix A 



4 An arbitrary basis of normal vectors is reconstructed from the normal basis {n, 1} by a local rotation: 

ni = cos £1 n + sin ill ; ri2 = — sin ft n + cos fl 1 . 



(11) 
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In the adapted frame, the normal connection u a u is identified as a surface vector, 



AB 



(15) 



where £ A are the components of the tangent vector e a with respect to the basis of tangent vectors to 
the hypersurface {E^}, e a = £ a Hja with £ A = dU A /ds a , and l A are the components of the vector 
normal to T N , tangent to Sjv+i, 1 = l A E A . 



Straightforward calculations give (details are presented in Appendix C) 



-V Z K + K [ — - K ab K ab + /i I + K {1) o 



(16a) 



(1) 



- K ab K ab{l) I +K, 



41) 



A (l) 



K a6(1) K o6 «+^ +fi n , (16b) 



where i^m a 6 is identified with the projection of Kab or| to T^at+i, and thus given by (see Appendix A) 

K(l)ab = ^a^b^AB ■ 

Vt\ and Vt n are given by 

+ \cu\ 2 n ; 

f2 n = —kW-uj — 2uj-'Wk+\ijj\ 2 K^i). 
They both vanish when u a does. The divergence and squared norm of cu a are given by 

V-w = (^/ B -^/ c V c )KAB-^ afc{1) ; 



(17) 



(18a) 
(18b) 



(19a) 
(19b) 

where H AB = g ab £ A £ B = G AB - l A l B is the projector onto rV,jv+i and R A b is the Ricci tensor on the 
hypersurface. 



2.1 Confinement of an elastic curve by a surface 



Before examining surfaces confined within a hypersphere, it is useful to confirm that Eqs.(16a) and (16b) 
reproduce the expressions describing the confinement of elastic curves by a two-dimensional surface in 
three-dimensional Euclidean space derived in Ref. [5]. 

Let the curve be parametrized by arc-length s. Its tangent e s = T is now a unit vector so that 
g ss = 1; V s is the derivative with respect to s, denoted by a prime '. The extrinsic curvature tensors 
for each of the two normals are scalars along the curve: one identifies the geodesic 

curvature; the projection Ksl = K^ = K n as the normal curvature; and the normal connection 
w s i2 = —t A l B Kab = T g as the geodesic torsion. On account of these identities, the normal covariant 
derivatives assume the form V S K^ = n' n + n g r g and V s k = n' g — n n T g . The stress tensor is now 
represented by a vector F along the curve, defined by 

F := f = ( K l±< _ ^ T _ ( K ' n + KgTg) n - (< - n n r g ) 1 . (20) 
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In addition, V • u — r' and |c<;| 2 = r 2 , so that Eqs.(16a) and (16b) reduce to 



- £i 



F' ■ 1 = K" + K g 



F' n = K n n + K n 



Kg + K 2 n 



y n T g ) 



i: 



(21a) 
(21b) 



Setting E\ = reproduces the Euler-Lagrange equation describing a confined elastic curve obtained by 



Nickerson and Manning [H|2]. Eq. (21b) identifies the magnitude of the normal force written down in 

If the curve is confined by a sphere (with K n = 1 and r g = 0), its equilibrium states are described by 
the Euler-Lagrange equation, —ey = K" g + K g (K 2 /2 — // + 1/2) [IB]. The corresponding normal force 
simplifies to A = — e n = k 2 /2 — /i + 1/2. A detailed treatment of confinement of elastic curves by 
spheres was presented in |6J. In this case, as shown in |6J, it is possible to exploit the residual rotational 
invariance of the problem. This implies that the torque vector M is conserved, where M is given by 
M = XxF + S, a sum of the moment of the force F and an intrinsic moment S = — K g n + K n lj^] The 
equation M 2 = constant provides a quadrature for K g . The loop is then constructed directly from its 
curvature data using a polar chart adapted to the vector M. 



3 Confinement by hyperspheres 

3.1 Hyperspherical shape equation and transmitted force 

If the hypersurface is a unit iV + 1-sphere, its extrinsic curvature is given by Kab = Gab so that 



K = N + 1 = K AB K AB - Using Eq.(17), one then identifies K t 



addition, u a = —E^Ia = 0, so that Oi = = Vt n . 
The stress tensor is given by f a = f ab e b + f a l, where 



L a6 



9ab 



and therefore K w = N. In 



f ab = K [ K ab - - K g ab 



N 



(N-2)+^)g ab ; 



The spherical Euler-Lagrange equation is then given by 

( K 2 N 

ei = V a f a • 1 = -V 2 k + k [ — - K ab K ab + —(N - 2) + fi 



The magnitude of the corresponding normal force is 



- A = s n = V a P ■ n 



k 2 + N 2 



(N - 2) + fj,N . 



(22a) 
(22b) 

(23) 



(24) 



These are the hyperspherical reductions of Eqs. (16a) and (16b). An alternative derivation of Eq.(23) is 
provided in |Appendix D which exploits the hyperspherical geometry in the variational principle. Notice, 



5 M is the one-dimensional reduction of the surface torque tensor, m al = e l ^ k : X 3 f ak : + S m , where S al = K 1 e a ^n 



discussed in Appendix E 
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however, that the derivation presented there suffers the disadvantage of having nothing to say about the 



normal forces described by Eq.(24). This involves the behavior of the surface energy under deformations 
normal to the hypersurface and thus falls beyond the scope of any derivation that is intrinsic to the 
hypersurface. 

When N = 2, the force A is completely determined by the multiplier /x, e n = 2/x. It is thus constant 
along the contact. This is not true if the confining surface is not spherical. In the next section, the 
origin of this coincidence will be traced to the scaling behavior of the constrained energy. 



3.2 Virial identities 



The bending energy of a two-dimensional surface has one particularly distinctive feature: its scale 
invariance. It is also invariant under the conformal transformations of the surfaces induced by embedding 
from its environment [7J. A consequence of scale invariance is that tension is not introduced in the surface 
when its area is fixed unless additional constraints are introduced to set a scale; a circular loop, on the 
other hand, would dilate if its length were not constrained; a three or higher dimensional surface would 
collapse. This behavior reflects the dimensional dependence of the behavior of bending energy under 
scaling. Confinement will break the scale invariance of the two-dimensional bending energy. 

More specifically, the behavior under scaling of the constrained energy in equilibrium implies the identity 

dA N n • Y A = — (iV — 2) H[Y] — fiNA N , (25) 



connecting a weighted surface averaged A, the bending energy H, and the multiplier \i. When N = 2, 
there is no explicit dependence on H. To derive Eq. (25), consider the effect of a rescaling Y — > AY 
on the constrained energy H c , given by Eq.((2]). This gives 

F C [AY, U\ = A N ~ 2 H[Y\ + A N J dA N X(s) ■ [AY(s) - X(U(s))] + fi(A N A N - Aq) , (26) 

so that in equilibrium, 



dH c [AY,U] 



dA 







A=l 



(27) 



This reproduces the identity Eq. (25). 

In the case of interest, N = 2, the bending energy is scale invariant; it does not appear explicitly in Eq. 

If the surface T N is closed, Eq. ([8]) implies the integrability condition 

dA N An = . 



(28) 

Contact with the hypersurface does not need to be complete. This identity also follows from the 
translation invariance of Eq. (25). More explicitly, let Y — > Y + a, where a is a constant vector. The 
only non-invariant term in Eq. (25) is the first one. 



If £jv +1 is a unit hypersphere, Eq. (25) implies the remarkable formula 

II := / dA 7v A = —(N — 2) H[Y] - fxNA 



N 



(29) 
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connecting the total force transmitted to the hypersphere to bending energy, the multiplier /i and the 



area. While it is also implied directly by Eq. (24), its derivation as an identity associated with scaling 
presented here does not require any local input, and as such is more compelling. There are two points 
worth emphasizing: 

(i) In the absence of a constraint on the area, there will be no force on the containing hypersphere: 
II is completely determined by the product fiA 2 . Indeed A itself is constant, proportional to /i. This 
is a direct consequence of the scale invariance of the bending energy of a two-dimensional surface. In 
striking contrast with the behavior of a confined circular loop, the confining force is due entirely to the 
area constraint and is constant over the surface. 



(ii) The contribution of bending energy to the transmitted force in Eq. (29) changes sign when N = 2. 
This reflects the contrasting behavior of a loop (iV = 1), and a three or higher dimensional surface. In 
the former case, the area (length) constraint works against the tendency of bending energy to expand 
the loop; whereas in the latter it prevents its collapse. 

3.3 Weakly Confined Spheres 

Consider a surface of area Aq/Ah — 1 <C 1 so that it is weakly confined within a unit hypersphere. In 
the linear approximation, the surface can be described in terms of a deformation of a two-sphere of unit 
radius (X say) along the normal 1 into the hypersphere. The surface can thus be characterized by the 
amplitude of the deformation 

X = X + $ 1 . (30) 

More generally, consider the deformation of a surface with metric tensor g ab and extrinsic curvature K ab 
embedded in a Riemannian manifold with metric Gab- The metric and extrinsic curvature of the surface 
deformed along the normal 1 are given by 5g ab = 2K a b$, and 

5K ab = -V a V 6 $ + [K ac K c b - R ABCD l A l c Ea^ , (31) 

where Rabcd is the Riemann tensor constructed with Gab (see, for example, [19J for a derivation). 
Thus 

5k = -V 2 $ - [K ab K ah + R A bI A 1 B }<$> ■ (32) 
The sphere is totally geodesically embedded in the hypersphere, so that n ab = 0. For a deformed unit 



sphere, Eq. (32) thus reads 

S«=-[Vg + 2]$, (33) 
where Vq is the Laplacian on the unit two-sphere. There is no corresponding first order change in the 



induced metric. The linearization of Eq. (23) is thus given by 

[V^ + 2 + /i][V^ + 2]$ = 0. (34) 
$ is now expanded in spherical harmonics, 



There are three zero modes with I = 1, which correspond to a rotated two-sphere. The mode with I = 
corresponds to a dilatation of the sphere changing its area so is excluded by the area constraint. Global 
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solutions of Eq.(34) imply a fixed single value of I. The linearized equilibrium states are thus spherical 
harmonic deformations of the two-sphere. This result is also directly analogous to the behavior of an 
elastic curve confined by a sphere |6J. The force per unit area transmitted to the hypersphere is given, 
for each I, by 

t i l = l(l + l)-2. (36) 

It is constant and independent of the magnitude of $. The latter observation may appear to be 
counter-intuitive: it can, however, be understood as a higher dimensional analogue of the Euler buckling 
instability in an unstretchable rod [9J. A given critical force will be necessary to buckle the sphere into 
each harmonic. The total force does, however, depend on the excess area. 

To determine the energy of the deformed sphere, note that the excess area is related to the deformation 
as follows 

A « An + j dQ 2 5k<5> = An - J dQ 2 $[Vg + 2]$ , (37) 
where cff2 2 is the area measure on the unit two-sphere. Thus 

A - An = fit \ a im\ 2 - (38) 

m 

The energy of the deformed state with a fixed value of I is then given by[[] 



Hi 



dtt 2 [5k 2 + A5k§] 
-(^-2)^ + 2)^2 K 



1 (^-2)^ + 2) 

2 a** 



(A - An) 



(39) 



where ^ is given by Eq(36). It increases linearly with area. Energy per unit excess area, Hi 
Hi/ (A — Air), is plotted as a function of / in Fig. [2] it increases monotonically with /. 




Figure 2: Energy per excess area of the deformed stated for low values of I 



6 The first order change in area, given by SA = [ d^nQ , vanishes [151 119], 

7 This is obtained most directly by expanding Eq. (D. 4| with TV = 2 about a unit two-sphere. 
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The ground states are represented by ellipsoidal (or quadrupole) deformations with I = 2. Modulo its 
orientation with respect to the sphere, the ground state is three fold degenerate]^] One would expect 
this degeneracy to be lifted at higher orders in the perturbation. 

Each of these perturbative states of the confined sphere will possess analogues in the full non-linear 
description. In particular, an infinite number of axially symmetric states solutions, one for each /, would 
be expected; ellipsoids, pears, and dumbbells (or discocytes) corresponding to / = 2,3,4 respectively. 
A detailed analysis is beyond the scope of this paper but will be taken up elsewhere [20] . It would be 
interesting to know if cross talk between different values of I emerges as non-perturbative features of 
the confinement process giving rise to composite structures such as (cup-shaped) stomatocytes. 

The analysis of the confined elastic loop suggests that the simple dependence of the energy and the 
transmitted forces on the excess area will diverge from the simple behavior described here as the confined 
surface explores more of the hyperspherical volume and deviations from geodesic behavior become 
increasing pronounced. 

In analogy with an elastic loop, it is reasonable to expect that the only stable weakly confined states 
are the ellipsoidal ground states with / = 2. All excited states (with I > 3) are unstable. To show this, 
note that the second variation of the energy about a weakly confined state is given by 



Ho 



(40) 



(41) 



where the operator Ci is defined by 

Cl = [Vl + l(l + l)}[Vl + 2]. 
The eigenvalues of Ci are identified as 

X k = (k(k + 1) - 1(1 + l))(k(k + 1) - 2) . (42) 
These eigenvalues are plotted for / = 2, 3, 4 and 5 in Fig. [3] One observes that for each / > 2 there will 

600 
400 
200 



-200 




Figure 3: Eigenvalues of C for low values of I. Each point corresponds to 2k + 1 deformation modes on account 
of the spherical harmonics degeneracy 

be I 2 — 4 negative eigenvalues corresponding to unstable modes of decay. Thus, as claimed previously, 
8 Note that, at this order, prolate and oblate deformations of the sphere are interchanged under a sign change in $. 
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the only stable states are those with I = 2. This generalizes the qualitatively very similar stability analysis 
of loops weakly confined within a sphere |6J. Whether this behavior is a faithful representation of the 
behavior when the surface is not weakly confined, on the other hand, is anyone's guess. Indeed, it will 
no longer be legitimate to assume that the ground state is completely attached. 

4 Conclusions 

In this paper, the confinement of a freely deformable closed two-dimensional surface within a three- 
dimensional hypersphere of smaller radius was examined. This is the natural generalization of a problem 
of a more domestic nature: the confinement of a loop within a surface. In the absence of the constraint 
the unique non-singular-topologically spherical-equilibrium state is a round sphere. Under confinement, 
this sphere deforms into one of an infinite number of equilibrium shapes; they have been examined in 
the limit of weak confinement. Mathematically, this problem can be viewed as a natural generalization 
of the spherical Willmore problem, to which it reduces when the constraint on the area is relaxed. In 
this context, it is also reasonable to expect that counterparts of these states exist in each topological 
sector of the two-dimensional surface. 

While the study of confinement in higher dimensions is not expected-or is unlikely— to offer insight into 
problems that arise in biology or soft matter, the questions posed are of a fundamental nature that have 
not been explored previously. Indeed, if the geometrical invariants of a surface are ordered by derivatives, 
after area, the bending energy is arguably the most important-non-topological-invariant among them. 
Analogues of the problem could be relevant in relativistic field theories, be it the Euclidean formulation 
of such theories, or the behavior of stationary spacelike surface states. Analogous actions were discussed 
by Polyakov as effective field theories [TU] , 

If the analysis of confined loops is any guide, a non-perturbative extension of this work should be worth 
pursuing. To this end, it is possible to exploit the conformal invariance of the two-dimensional bending 
energy |7] to reformulate the problem in a rather striking way. For under inversion in a hypersphere, 
itself located on the confining hypersphere, the latter will get mapped to a three-dimensional Euclidean 
hyperplane. The confined surface on the three-sphere is replaced by a surface that is free to bend in 
this three-dimensional space without obstruction. The area, however, is not invariant under conformal 
transformations; it gets replaced by a surface integral 



involving a potential that will depend on the Euclidean distance |X| from points on the surface to the 
origin of the three dimensional space. The confining surface is replaced by this attractive potential. 

It is also possible to consider confinement as a dynamical process in a relativistic theory of extended 
objects. Two-dimensional surfaces are interpreted now as the world-sheets of strings |2T] 122]; the 
confining hypersphere is replaced by a three-dimensional de Sitter space, the fixed worldsheet of a 
membrane. 




(43) 
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Appendix A Embedding identities 

In this appendix, several useful identities connecting Tn,Tn,n+i and T, N+1 , as defined in section |2j are 
collected. 

Let Eia, A = 1, . . . N + 1 be the tangent vectors to the hypersurface Ejv+i, and n its normal vector. 
The Gauss structure equations describing the embedding of S^+i in Euclidean space are given by 



V/iE 



-K AB n. . 



(A.l) 



Here is the covariant derivative on Sat+i compatible with its induced metric, Gab = E^ ■ E B . 
These equations define the curvature tensor Kab- 



The counterpart of Eq.( A.l ), describing the embedding of the surface Tn directly into the same Euclidean 
space, is given by 

V a e b = -K abI n I . (A.2) 

Choose ri! = n, and n 2 = 1. The surface tangent vectors e a , a — 1, ■ ■ ■ , N can be cast as a linear 
combination of their hypersurface counterparts as follows, e a = £ a ~Ea, where 8 A = dU A /ds a . The 
vector 1 can likewise be expanded, 1 = I a Ea- 

Projecting V^E^ onto gives 



£ a £ b V a Vb = V a e b - (V«#) E B . 



(A.3) 



It follows from the structure equation (A.l) and its definition, that the extrinsic curvature of T N asso- 
ciated with the rotations of n onto is given by Eq.( 17 ), K ab (i) = £ A £ B Kab- In particular, using 
the completeness of the basis vectors {£ A ,l A }, g ab £ A £f + l A l B = G AB , the surface trace of Kn\ ab is 
given by 

K {1) = g ah E A E B K ab = (G AB - l A l B )K AB =K- l A l B K AB • (A.4) 

In addition, 

E A l B V A E B = V a l - (VJ B )E B , (A.5) 

so that 

0J a i2 ■= n • VJ = -E A l B K AB ■ (A.6) 



If the surface is embedded in a hypersphere with Kab = Gab, then K^ ab = g ab and cu al 2 = 0. 

For a curve, with N = 1, with unit tangent vector T = t A F,A, Km = t A t B KAB = K n ; ^12 
l A t B KAB = T g , where n n and r g are the normal curvature and geodesic torsion of the curve. 



The counterpart of Eqs.(A.l) and (A.2) for T 



N,N+1 



embedded into T, N+ i are 



D f A 



(A.7) 
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where 

D a S° := V £ 6 C + T C AB £X B , (A.8) 

and T^b ' s tne Christoffel connection compatible with the induced metric on Sjv+i, Gab- Eq.(A.7) 
defines the extrinsic curvature tensor n ab associated with the rotation of the normal 1 onto Tn. This 

(2) 

provides the identity for K ab , 

:=K ab = -l B D a £* . (A.9) 



Appendix B Derivation of the stress tensor 

In this appendix, the stress tensor introduced in Eq.(|3]) will be derived. This will be done by extending 
to co-dimensions higher than one the auxiliary variables method introduced in [IT]. One constructs 
a functional H C [X., e a , n, g ab , K ab , co a IJ , f a , X a j, A/j, X ab , A ab T , A a u ] treating the embedding functions 
Y, the adapted basis {e^n 7 }, the metric tensor g ab , the extrinsic curvature tensor K* b , as well as the 
connection u a IJ as independent variables, by introducing Lagrange multipliers f a , A a J; Ajj, X ab , A ab j, A a 7J 




f° ■ (e„ - d a Y) + X ab (g ab - e a ■ e b ) + ^X U (5 IJ - n 1 ■ n J ) 



H C = ^J dA N Ki K 1 + J dAj 

+ J dA N [X%e a ■ n 1 + A^j^Kj - d a n J ■ e b ) + A a jj(u; a IJ - n 1 ■ d a ii J )] . (B. 1) 

The bending energy now depends only on the independent variables g ab and K* b . Lagrange multipliers 
with more than one index possess the same symmetry as the quantities they multiply. Thus, for instance, 
X ab = X ba , Xu = Xji, A a jj = -A a jj. The variation of H with respect to Y gives 

5*H C = [ dA N VJ a -5Y. (B. 2) 



Thus, in equilibrium, f" will be conserved on a free surface; it is identified as the stress tensor on the 
surface [17J. The Euler-Lagrange equations for e a provide and expansion of the Lagrange multiplier f a 
with respect to the adapted basis 

r = (2X ab + A acI K c b j) e b - (X aI - A abJ u b /) n T . (B. 3) 

The Euler-Lagrange equation for n is 

(A a/ + V b A abI - 2A bIJ K b a j)e a - (V a A a/J + A abI K ab J + A /J )nj = . (B. 4) 

Thus from the linear independence of the adapted basis, one identifies 

\ aI = -V fe A a6/ + 2A bIJ K b a j , X IJ = -V a A aIJ - A abI K ab J . (B. 5) 

There remains to determine the three Lagrange multipliers X ab , A abI and A aIJ . They are determined 
from the Euler-Lagrange equations for g ab , K ab j and co a IJ respectively: 

X ab = K I (K abI - ^AfV b ) ; A abI = -AT 7 ^ ab ; A aIJ = . (B. 6) 



9 The Christoffel connection r c ab is constructed using g a b and its derivatives; in contrast, the normal connection Lu a IJ 
is not constructed from the normal "metric" 5 IJ so that, in general, it is necessary to implement its definition in the 
variational principle. 
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A aIJ vanishes because the connection cu a IJ does not appear explicitly in the bending energy. Had the 
definition of ui a IJ been omitted, no error would have been incurred. However, one would not be so lucky 
if the energy had depended explicitly on co a IJ . 



Substituting the expressions for the Lagrange multipliers given by Eqs. (B. 5) and (B. 6) into Eq.(B 
3|, f a assumes the form 

1 



r = K I \K 



-abl 



) e 6 - VTn; 



a. rsl. 



(B. 7) 



where V a is the covariant derivative involving the normal connection tujj 



Appendix C Hypersurface adapted Euler-Lagrange deriva- 
tives 



In this appendix, the details of the derivation of Eqs. (16a) and (16b) are provided. First dismantle the 
covariant Laplacian V 2 K ! in terms of the normal connection, so that 



W 2 K T 



jib 



9 I » a 

(VbKl) — UaljVbKj 

= V 2 Kj + V a u aIJ K J + 2oo aIJ V a K J + u aI J u aJL K L 
The components along n and 1 are given respectively by 

V 2 K (1) = V 2 K {1) + kV ■ u + 2uo ■ Vk - \uj\ 2 K ( 



^(i) 



(C. 1) 

(C. 2a) 
(C. 2b) 

where V • u — V a u) a i2 is the divergence of the normal connection and \co\ 2 = u a 12 co a 12 is its squared 
modulus. For an embedding in a hypersphere, the normal connection vanishes, cu a i 2 = 0, so that the 
normal covariant derivatives are replaced by covariant derivatives, V 2 k = V 2 k; in addition, Km is 
constant, so that V 2 K^ = 0. 



V 2 K 



V 2 k - Km V • u - 2uj ■ VKm - \oj\ 2 k . 



More generally, Eq. (A. 6) is used to express the divergence and squared modulus of the normal connection 



in terms of hypersurface curvatures. The second term is given by 

V-w: = V a u n ™ = -V a (££l B K AB ) 



■>al2 

Kl A l B K AB - K ab £*£ b B K AB - g ac £^l B V c K AB 

H AB l c V c K AB . 



k\ a \ b K ab 



K ab K, 



ab(l) 



(C. 3) 



Here 1-L AB = g ab £ A £ b B represents the projector onto IVat+i- On the second line, the relations V a l B = 
K ab £ B and D a £ A = —nl A are used, which follow from the structure equations associated with the 
embedding of the surface IVtv+i i n ^n+i< (Eq- (A. 9) and its contraction). In addition, the Codazzi- 



Mainardi equation for Ejv+i, V ' A K B c = V B K A c has been used in the third line. 



The squared modulus is given by 

M 2 : 



9 ab oo al2 io bl 2 = g ab £ A £ b B l c l D K AC K BD 
l A l B [K AC K B C - K AC K BD l c l D ) 
l A l B [K(i)K AB — Rab) i 



(C. 4) 
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where the completeness of the basis {£ A ,l A } and the contracted Gauss-Codazzi equation have been 
used. Here Rab is the Ricci tensor constructed with the metric Gab- 

One also requires the decomposition of Kjf ab where f ab = f a ■ e b = Ki{K abl - l/2g ab K I ), which are 
given by 



(1) rab 



K K 'f 

^ab J 



U ab{x) K^ - l -K 2 {1) -»\ + k L ab K^ - \kK<A ; (C. 5a) 

K ab f ab = L ab K ab W - l -KK^ + K L ab K ab - l -K 2 - ^ . (C. 5b) 



(1) 



For a surface embedded in a hypersphere, with K t 

and (K ab{1) K ab ^ - K 2 {1) /2 \ = -N 2 (N - 2)/2. 



g ab , k {n ab K^ ab - kK^/2) = -(N-2)k 2 /2 



For confined curves, with = n n and k = K g , Eq, (C. 5a) reduces to K n [(K 2 + «£)/2 — fj], whereas 
Eq.(C. 5b) reduces to K g [(K 2 + k 2 )/2 — fi]. 



Appendix D Variational principle adapted to a hypersphere 

It is instructive to decompose the bending energy into parts adapted to its confined environment. Using 



Eqs.flA^D and <\A^ one finds 

H = \ JdA N K^K {I) = X - JdA N [k 2 + (K- l A l B K AB ) 2 } ■ (D. 1) 

The first term is the bending energy of the surface embedded into the three-dimensional curved geometry 
described by the metric tensor Gab, 

Hv= l - j dA N n 2 . (D. 2) 

The second is the bending energy inherited from the hypersurface. This is the generalization to higher 
dimensions of the Pythagorean decomposition of the Frenet curvature k of a space curve into geodesic 
and normal parts, k 2 = k 2 + familiar in the theory of surfaces [5]. 



The Euler Lagrange equation for the surface s — > U(s) described by the bending energy (D. 2) is given 
by 

- V 2 k + - (k 2 - 2n ab K ab )K - R AB l A l B n + hk = . (D. 3) 

This equation was derived in the 90s in the context of relativistic extended objects [19J. 

For an iV + 1-sphere, the Ricci tensor is proportional to the metric tensor, Rab = NGab 50 that 
RabI A 1 B = In this case K A b = G A b, and K = N + 1 so that Eq. flD. l[ ) then reduces to 



H = i J dA N [k 2 + N 2 



(D.4) 



10 The Gauss-Codazzi equations describing the embedding of the hyperface in the Euclidean background, Rab — 
KKab — KacK b, together with the identity, K ab = Gab, are useful for remembering the relationship between Rab 
and Gab- Of course this relationship is completely intrinsic. 
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The normal curvature contributes a constant energy. If the Euler-Lagrange derivative corresponding to 
this contribution (N 2 k/2) is added on the left hand side of Eq. (D. 3), the shape equation Eq. (23) is 
reproduced. 



Appendix E Rotational invariance and confinement by a hy- 

persphere 

If the confining surface is a hypersphere, the Euler Lagrange equations will respect rotational invariance 
in IR 4 . Thus, whereas stress is not conserved, torques will be, corresponding as they do to the Noether 
current associated with the rotational symmetry of the energy. 

Consider a small rotation about the origin 5X l = e yW fi jfe Jf', characterized by the antisymmetric rotation 
matrix f2 y , where e^ kl is the Levi-Civita symbol and k, I — 1, . . . , 4 index tensors in Euclidean space 
M 4 . Each of the two normal vectors n 7 ,1 = 1,2 rotates accordingly: 5n l1 = e^ kl ^V k n 11 . It can be 
shown that, in equilibrium, on any surface patch 

5H C = J dA 2 V a m ija , (E. 1) 

where 

m aij = e ijkl ^ak x l _ Kig ab e k n Ilj _ (E. 2) 

Here f ak = f a • x fc , e k = e a ■ x fc , etc. 

Thus m l, ' a is conserved if H c is rotationally invariant. This is confirmed by the following calculation: 

V a m aij = e ijkl (v a f ak X l + f ak e l a - K^^e^ - V a K ig ah e k n 11 - K^jn^n 11 ^ . (E. 3) 

In general, all but the first term either vanish identically or cancel among themselves. For the first term 
one has e l ^ kl V a f ak X l = e^ kl e a n k X l . However, if the surface is confined to a unit hypersphere centered 
on the origin, so that X = n, this term also vanishes. 

It is not obvious how to exploit this invariance in a manner analogous to that in the treatment of a 
polymer confined within a sphere. In any case a better strategy is to exploit the conformal invariance of 
the two-dimensional bending energy, as sketched in section |4j 
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